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摘 要 :在 输 水 工程 中 ,由 于 地 形 、 地 质 等 因素 的 限制 ,分 又 型 结构 管道 得 以 广泛 推广 和 应 用 ,分 又 
结构 输 水 管道 泄漏 是 国内 外 学 者 非常 关注 的 一 个 重要 问题 ,而 瞬 交 流 法 是 目前 高 效 、 精 确 且 快速 的 
泄漏 检测 方法 之 一 。 基 于 此 ,本 研究 构建 了 含 分 又 结构 输 水 管道 的 肯 变 流 模 型 ,探究 了 影响 阀门 末 
端 压 力 变化 的 因素 。 结 果 表 明 ,压力 曲线 的 差异 性 随 分 叉 管 长 度 的 增加 而 增 大 ;压力 曲线 的 衰减 速 
率 随 阀门 关闭 时 间 的 增加 而 变 慢 ;泄漏 量 对 压力 曲线 变化 周期 影响 不 大 ,但 随 着 泄漏 量 增 大 ,压力 
曲线 的 面积 积分 变 小 ;泄漏 位 置 不 同 ,压力 曲线 在 幅 值 处 的 形状 不 同 , 当 泄 漏 位 置 位 于 管道 上 半 段 
时 ,曲线 在 幅 值 处 会 凸 起 , 当 位 于 管道 下 半 段 时 ,曲线 在 幅 值 处 会 出 现 凹 陷 的 趋势 。 
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Numerical simulation of leakage transient characteristics 
of water pipeline containing bifurcation structure 
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Abstract :Im water transportation projects , due to the constraints of topography and geology, bifurcated 
pipelines have been widely promoted and applied. The leakage of bifurcated water pipelines is an important 
jssue of great concern to scholars at home and abroad ,and the transient flow method is one of the most effi- 
cient ,accurate and fast leak detection methods at present. Based on this ,this paper constructs a transient 
fow model of water pipeline with bifurcation structure ,and explores the influence of different factors on the 
pressure change at the end of the valve. The factors affecting the pressure change at the end of the valve 
were investigated. The results show that the difference of pressure curve increases with the increase of the 
length of bifurcation pipe ,the attenuation rate of pressure curve decreases with the increase of valve clos- 
ing time ,and the leakage has little influence on the variation period of pressure curve. The leakage position 


is different ,the shape of the pressure curve at the amplitude is different. When the leakage position is in 
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the upper half of the pipeline ,the curve will be convex at the amplitude ,and when it is in the lower half of 


the pipeline ,the curve will be concave at the amplitude trend. 


Key words :transient flow; numerical simulation ;bifurcated pipeline ;leak detection 


随 着 社会 经 济 的 快速 发 展 ,人 们 对 水 资源 的 需 
求 不 断 增加 ,由 于 长 距离 输 水 系统 常常 经 过 复杂 地 
形 且 支 路 多 ,实际 工程 中 广泛 应 用 分 又 结构 管道 。 
然而 ,由 于 腐蚀 老化 .人 为 破坏 等 原因 而 导致 管道 泄 
漏 问 题 越 来 越 严重 ,威胁 了 输 水 系统 的 安全 运行 "1。 
管道 泄漏 不 仅 会 造成 严重 的 经 济 损失 ,也 在 一 定 程 
度 上 对 环境 造成 污染 ” ,甚至 会 对 人 民 的 生命 健康 
造成 威胁 ,因此 ,迫切 需要 提出 一 种 精确 高 效 的 管道 
泄漏 检测 方法 。 瞬 变 流 法 是 检测 管道 泄漏 的 有 效 手 
段 之 一 中 ,其 主要 原理 是 :在 管内 施加 一 个 激励 信号 
来 产生 水 击 波 , 当 该 水 击 波 经 过 泄漏 点 时 会 发 生 误 
减 和 畸变 ,由 其 变化 规律 及 特征 可 以 判断 泄漏 位 置 
hh 浊 漏 量 。 
一 目前 ,针对 输 水 管道 泄漏 问题 ,国内 外 学 者 进行 
也 系列 的 研究 。 陈 玉 涛 等 中 构建 了 一 维 数学 模 
到 办 析 了 影响 分 又 管道 压力 变化 的 因素 。 周 知 进 
等 党 研究 了 不 同 角度 下 流体 冲击 对 管 壁 应 力 的 影 
嘛 N 蒋 明 等 "发 现 应 用 特征 法 求解 克 合 瞬 变 流 问 题 
会 造成 较 大 的 误差 ,提出 了 一 种 求解 管道 耦合 瞬 变 
模型 的 计算 格式 ,并 通过 实验 与 仿真 证 实 该 方法 的 
可 每 性 与 有 效 性 。 赵 点 等 “采用 断 电 停 泵 方式 引入 
瞬 需 流 ,分 析 空 气 阀 与 空气 负 结 合 对 水 锤 的 防护 效 
果 6 发 现 其 大 大 降低 了 管道 的 水 锤 压力 。 陈 兵 等 
利 甫 0LCA 软件 模拟 了 泄漏 过 程 中 压力 变化 特征 ， 
并 综 析 了 影响 因素 ,为 管道 输送 的 安全 设计 与 运行 
提出 了 合理 建议 。 白 莉 等 ”基于 瞬 变 流 理论 ,用 非 
线性 动态 系统 的 状态 空间 模型 来 蔡 代 水 动力 学 模 
型 ,用 滤波 来 估计 预测 的 水 力 要 素 ,泄漏 定位 的 精度 
得 以 提高 。 蔡 正 敏 等 (" 结合 序 贯 概率 比 检验 法 与 
管道 泄漏 故障 诊断 ,提出 了 沪 漏 识别 因子 这 一 概念 ， 
管道 泄漏 检测 的 精度 得 以 提升 ,引入 了 瞬 变 流动 ,并 
建立 了 瞬 变 流 方程 ,对 不 同 泄漏 位 置 处 的 瞬 变 特性 
进行 了 对 比分 析 。 杨 开 林 、 郭 新 一 等 学 者 "中 针对 
长 距离 输 水 管线 构建 了 水 击 频 域 数学 模型 ,利用 快 
速 关闭 阀门 以 及 瞬 变 流 发 生 器 的 方式 引入 瞬 变 流 ， 
泄漏 位 置 由 瞬 变 信号 的 尖 点 来 确定 ,并 利用 遗传 算 
法 求解 了 反问 题 ,对 实际 工程 中 泄漏 检测 的 精准 定 
位 提供 了 可 能 。 杨 镇 榜 等 "提出 了 一 种 基于 瞬 变 
流 理论 的 输油管 道 泄漏 检测 方法 ,以 连续 压力 波 为 


输入 信号 ,以 管道 出 口 压力 信号 为 输出 信号 。 李 六 
佳 等 "通过 瞬 变 流 泄 漏 仿 真 模型 来 模拟 压力 波形 ， 
并 建成 了 模式 识别 库 ,管道 是 否 发 生 泄漏 由 负 压 波 
识别 法 来 判断 。 陈 法 等 '" 建立 了 输 气管 道 瞬 变 流 
模型 ,分 析 了 发 生 泄漏 时 管道 沿 程 置 压力 ,流量 随时 
间 的 变化 特征 ,并 确定 了 泄漏 位 置 。 输 水 过 程 中 由 
于 管道 泄漏 ,阀门 关闭 等 引起 瞬 态 流动 时 ,流动 信号 
在 噪声 ,振动 等 的 干扰 下 ,容易 引发 错误 报警 ,因此 ， 
郑 杰 等 “采用 信号 采集 方法 设计 了 一 种 泄漏 检测 
装置 ,为 精确 检测 瞬 变 特性 提供 了 一 定 的 理论 指导 
意义 。 康 完 芝 等 2 采用 Flowmaster 软件 对 长 距离 
输 水 管道 进行 了 数值 仿真 计算 , 谐 波 分 量 根据 压力 
言 号 的 时 频 分 析 来 提取 ,并 根据 其 训 减 值 定位 泄漏 
位 置 。 孙 良 等 "提出 了 一 种 采用 特征 线 法 进行 汇 
漏 检 测 与 定位 仿真 模拟 的 瞬 变 流 模型 。 胡 靖 '”” 通 
过 用 反问 题 来 分 析 瞬 变 流 模型 的 方式 ,研究 分 析 了 
长 距离 输 水 管道 的 泄漏 。 

DUAN 等 .综述 性 介绍 了 目前 瞬 变 流 理论 的 
研究 进展 ,并 合理 建议 了 该 领域 的 未 来 主要 工作 。 
吴 猛 猛 等 基于 瞬 变 流 理论 ,建立 了 典型 的 管道 汇 
漏 系统 模型 ,模拟 分 析 了 泄漏 发 生 时 管道 瞬 态 压力 
波 的 影响 因素 。HE 等 ' 引 根据 瞬 变 压力 的 大 小 , 建 
立 了 3 种 基于 瞬 变 流 压 力 大 小 的 泄漏 量 计算 模型 ， 
分 析 了 影响 泄漏 量 值 的 临界 参数 ,并 通过 实验 验证 
了 模型 正确 性 。CHEN 等 构建 了 一 种 计算 管道 
瞬 态 泄漏 量 模型 ,求解 得 到 最 小 可 检测 泄漏 孔 的 位 
置 , 并 详细 分 析 了 泄漏 量 的 敏感 性 。WANG 等 ”1 通 
过 多 次 实验 发 现 噪声 影响 管道 泄漏 检测 的 精度 , 提 
出 了 滤 除 噪声 信号 干扰 的 方法 。 俞 逸 飞 等 ”通过 
数值 仿真 软件 Flowmaster 对 管道 泄漏 位 置 、 漆 漏 点 
大 小 、 直 径 变 化 等 进行 了 仿真 ,并 分 析 了 不 同 因素 对 
测量 精度 的 影响 。LEE 等 “建立 了 “ 直 管 -水 库 ” 
模型 ,通过 模型 试验 研究 了 因 阀 门 关闭 引发 的 瞬 变 
流 特 性 ,并 通过 逆 瞬 态 分 析 法 确定 了 泄漏 发 生 的 位 
置 以 及 泄漏 量 。 

至 今 ,已 有 一 些 国 内 外 学 者 对 管道 泄漏 时 的 瞬 
态 特 性 进行 了 研究 ,但 是 关于 复杂 分 又 输 水 管道 泄 
漏 的 研究 还 较 少 。 基 于 此 ,本 研究 针对 分 又 结构 输 
水 管道 ,构建 了 瞬 变 流 模型 ,模拟 分 析 了 其 泄漏 过 程 
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中 的 瞬 态 特性 变化 特征 ,分 析 Je 
闭 时 间 ` 汽 漏 量 泄漏 位 置 等 因素 对 管内 流体 瞬 变 
程 的 影响 规律 ,人 研究 成 果 可 为 长 距离 输 水 管 ne 
全 运行 提供 一 定 的 理论 指导 。 


1 计算 模型 


1.1 模型 建立 


分 又 结构 输 水 管道 泄漏 示意 图 如 图 1 所 示 , 由 
上 游 水 库 、 直 管 .分 叉 管 .下 游 水库 、 阀 门 等 组 成 。 假 

定 管道 为 一 均 质 直 管 ,管道 平 铺 , 不 考虑 倾斜 角 a 的 
影响 ,上 游 水 库 的 水 位 瓦 =25 m, 保 持 恒 定 不 变 , 上 
人 
;= 500 m, 管 道 直径 D =0.2 m, 流量 0。 = 
ms, 网 格 特征 长 度 Ax = 10 m, 时 间 步 长 Ar = 
es a 和 选择 ， 


上 游 水 库 
下 游 水 库 1 
下 游 阀 门 2 
上 游 直 管 天- 二- 


下 游 水 库 2 
图 1 含 分 义 结构 管道 泄漏 示意 图 


Fig. 1 


Schematic diagram of pipeline leakage 


with bifurcated structure 


1.23 控制 方程 


联 立 连续 方程 和 运动 方程 ,构建 了 输 水 管道 
变 流 数学 模型 , 即 
oH oH a 
万 = 一 二 1 =0， 
dt Ox go 
oH 1 vo folv| ov Ov 
mn gar 257 Be | 
(1) 
其 中 :k 为 Brunone 阻力 系数 ;f 为 达 西 摩 阻 系数 ;vw 


为 管内 流体 平均 流速 ;a 为 水 击 波 速 ;a 为 管道 倾斜 
角 ; 瑟 为 压力 ;g 为 重力 加 速度 ;x% 为 管道 位 置 ;i 为 
时 间 。 
1.3 ”特征 线 方程 

由 于 流体 流速 v 远 远 


元 小 于 水 击 波束 a, 因 此 模型 


中 "与 互 的 对 流 项 忽略 不 计 , 采 用 特征 线 法 将 式 
(1) 转 换 为 常 微分 方程 ,得 到 了 正 负 特征 线 方程 , 即 


+ a fa OO 
C a 
dx_ a 
dt 1+k 
-pa(l+k) 应 OO 
C- :dH 了 dQ De 
d 
下 号 (2) 
Dt 
将 B= 译 ,R= 5 ,010|= = Ooloul， 


Co = 1 + , 带 入 式 (2) ,并 沿 着 图 2 中 的 两 条 特征 线 
进行 积分 ,得 到 的 有 限 差 分 方程 如 下 
C*:H» - Hr + B(Qp - Qn) + R/CoQp | OQ. | = 0， 
C :Hs -Hs -BCo( Ob Q8) — ROb, |O, | =0 
(3) 
特征 线 网 格 节 点 处 示意 图 如 图 2 所 示 。 


图 2 节点 处 特征 线 网 格 
Fig.2 Feature line mesh at node 
将 整个 输 水 管道 平均 分 成 为 Y 段 , 每 段 的 长 度 
为 Ax ,时 间 步 长 为 At, 可 得 到 特征 线 插值 网 格 图 ,如 
图 3 所 示 。 


Ax A A 


图 3 管道 特征 线 插值 方式 网 格 图 
Fig.3 Grid diagram of pipeline characteristic line interpolation 
由 4 和 C 两 点 线性 插值 来 得 到 RR 点 的 压力 Hi 

和 流量 Qi ,如 式 (4) 所 示 。 
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Hr = (H, +KT DZCo， 

Qn = (0 + kQe) /Co (4) 
将 41 = Hi + BQi ,4, = B+R/C,|0,|, 
Hs - BCoQs ,D, = BCo +R|Q,| ,代入 式 (3)， 


司马 
纹 1 


CPP =4-40n， 
C-:H, = D, + D,0,» (5) 
联 立 求 解 式 (5) 中 C 和 C ,可 得 到 无 泄漏 管 
道 实时 的 压力 和 流量 ,计算 式 为 


_4 —D 
Cr A, +D,’ 
(4 + D, +D, -A,) 
H, = 1 本 2 0, (6) 


1. 参 边界 条 件 


1CD1 上 游 边界 条 件 
保持 上 游 水 位 恒定 , 即 万 为 定 值 ,假设 H=h,， 
代 六 负 特征 线 方程 中 可 得 上 游 节点 处 的 流量 为 
H,-D 
J (7) 


154232 泄漏 孔 处 边界 条 件 

也 假设 泄漏 孔 外 为 大 气压 ,泄漏 孔 处 的 泄漏 流量 
汶 记 ,泄漏 孔 上 游 流 量 为 0 ,泄漏 孔 下 游 流量 》 
0,>, 汇 漏 孔 处 的 泄漏 系数 为 C, ,泄漏 孔 处 边界 条 件 
的 示意 图 如 图 4 所 示 。 


图 4 ”泄漏 孔 处 的 边界 条 件 
Fig.4 Boundary conditions at leakage holes 
根据 质量 守恒 定律 ,可 得 
Q7 = Qi + (8) 
假定 泄漏 孔 处 流体 的 出 流 方式 为 小 孔 出 流 , 则 
可 求 得 泄漏 孔 处 的 泄漏 流量 0, 为 


Qi = C。 V2gH; (9) 
联 立正 、 负 特征 线 方程 可 计算 出 泄漏 点 处 压力 
及 上 、 下 游 流量 的 表达 式 为 
ee D,A! + A,D, - A,D,D; 
. D, + 4， 
A H, -DD 
D, (10) 


Qi = 


一 万 
A, *, 0 


> 
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1.4.3 下 游 阀 门 处 边界 条 件 
假定 阀门 处 为 孔 口 出 流 形式 ,阀门 处 的 流量 计 


A, (Onr)” A, (Qn7) 了 4 (Onr) 
ws | 5 


2H, 


(11) 
其 中 :i 为 定常 流动 时 阀门 全 开 时 阀门 处 的 压力 ; 
Qi 为 定常 流动 时 阀门 全 开 时 闪 门 处 的 流量 ; 7 为 阀 
门 相对 开 度 。 
把 式 (10) 代 入 正 特征 线 方程 中 ,得 到 阀门 未 端 

的 压力 为 


Hpy = A, — A,Qpy 
1.4.4 分 又 管 处 的 边界 条 件 
分 义 管 处 的 节点 示意 图 ,如 图 5 所 示 。 


(12) 


下 游 分 又 2 


图 5 分 义 点 处 节点 示意 图 
Fig.5$ Schematic diagram of nodes at bifurcation points 
忽略 分 义 管 连接 处 的 能 量 损 失 ,假设 该 处 有 相 
同 的 水 头 , 上 游 直 管 采用 C* 特征 线 相 容 方程 ,下 游 
两 分 义 管 均 采用 C -特征 线 相 容 方程 ,得 到 分 又 管 节 
点 的 压力 为 


Ha;’ 二 Hb’ 三 He = (13) 


Aa; ' Do, * De, 
然后 再 将 压力 瑟 代 入 到 各 个 相 容 性 方程 , 反 推 
可 得 到 各 个 管道 节点 处 的 流量 @。 


2 模型 验证 


目前 ,关于 分 又 型 输 水 管道 泄漏 的 实验 数据 较 
少 ,但 分 又 管 的 控制 方程 与 直 管 控制 方程 相同 ,区 别 
仅 在 于 连接 点 处 的 流量 变化 。 因 此 ,本 研究 将 直 管 
泄漏 数值 计算 数据 与 郭 新 董 等 ” 的 实验 数据 进行 
了 对 比 , 结 果 如 图 6(a) 所 示 。 同 时 ,为 了 进一步 提 
高 数值 模型 的 可 靠 性 ,将 Flowmaster 软件 模拟 结 
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与 模型 计算 结果 进行 了 对 比 ,如 图 6(b) 所 示 。 

实验 模型 的 基本 参数 为 :初始 流量 Q, = 10% 
Qo; 非 恒定 摩 阻 系数 为 0.014; 汇 漏 位 置 X; =750 m; 
其 他 参数 均 保持 不 变 。 
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(a) 直 管 阀门 末端 压力 变化 


一 一 编程 计算 结果 
-一 -Flowmaster 模 拟 结果 
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(b) 分 叉 管 阀门 末端 压力 变化 


图 6 验证 结 

Fig.6 Verification result 

"== 由 图 6 可 以 看 出 ,模型 计算 结果 与 实验 结果 符 
合 较 好 ,压力 衰减 趋势 及 变化 规律 保持 较 高 的 一 致 
性 。 两 者 之 间 的 平均 误差 约 为 7% ,最 大 误差 约 为 
8.4% ,主要 发 生 在 压力 峰值 处 ,这 是 由 于 不 同 模型 
计算 处 理 过 程 中 存在 差异 ,但 误差 均 在 允许 的 范围 
之 内 。 


naXiv 


3 计算 结果 与 分 析 


3.1 分 又 管 长 度 的 影响 


本 研究 分 析 了 3 种 不 同 长 度 分 又 管 对 泄漏 时 管 
内 瞬 变 流 特性 的 影响 ,长 度 分 别 为 250 .500 1 000 m， 
泄漏 量 为 0, 的 1% ,泄漏 位 置 位 于 距 上 游 水 库 500 m 
处 ,保持 其 他 参数 不 变 。 计 算得 到 了 不 同 分 又 管 长 
度 下 阀门 处 压力 的 变化 曲线 ,如 图 7 所 示 。 
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一 e 一 下 游 分 又 管 250 m 
一 一 下 游 分 又 管 500 m 
一 一 下 游 分 又 管 1 000 m 


ts 
图 7 分 义 管 长 度 对 下 游 阀 门 未 端 压力 变化 的 影响 


Fig.7 Pressure changes at the end of downstream 


valves with bifurcation length 

由 图 7 可 知 ,分 又 管 长 度 增加 2 倍 时 ,压力 峰值 
并 没有 出 现 完全 线性 变化 ,第 一 个 周期 处 的 压力 峰 
值 大 小 分 别 为 114. 86 .110. 97 .104. 78 m ,数值 相差 
越 来 越 大 ,周期 也 随 着 分 又 管 长 度 的 增加 而 变 长 。 
随 着 分 又 管 长 度 的 增加 ,压力 峰值 逐渐 降低 ,并 且 趋 
于 平缓 ,衰减 速率 逐渐 减 小 ,但 压力 曲线 的 变化 趋势 
基本 一 致 ,因此 可 知 不 同 分 叉 管 长 度 对 阀门 末端 压 
力 峰 值 的 影响 较 大 。 


3.2 阀门 关闭 时 间 的 影响 


研究 了 不 同 阀 门 关闭 时 间 对 阀门 末端 压力 变化 
规律 的 影响 ,下 游 分 又 管 长 度 为 500 m, 保 持 其 他 参 
数 不 变 ,阀门 的 关闭 时 间 为 0.5、2、5 s, 关 闭 方式 采 
取 线 性 关闭 的 方式 ,阀门 末端 压力 变化 曲线 如 图 8 
所 示 。 
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到 8 阀门 关闭 时 间 对 分 又 管 末端 压力 变化 的 影响 


Fig.8 The effect of valve closing time on the pressure 


change at the end of the bifurcated pipe 
由 图 8 可 知 , 随 着 阀门 关闭 时 间 的 增加 ,压力 曲 
线 峰 值 逐 渐 减 小 , 尖 点 逐渐 变 圆 滑 ,周期 越 长 ,压力 
曲线 差别 越 大 ,衰减 速率 也 越 慢 ; 当 阀门 关闭 时 间 为 
0.5.2s 时 ,达到 水 击 峰 值 的 时 间 为 3 s, 周 期 分 别 为 
7.1、8.0 s, 关 闭 时 间 为 5s 时 ,达到 峰值 的 时 间 为 5 s， 
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周期 为 8.5 s, 该 工 况 下 所 产生 的 水 击 都 是 直接 水 
击 , 各 压力 峰值 间 的 最 大 差 值 为 23.23 mm, 并 且 差 值 
随时 间 而 逐渐 减 小 ,说 明 阀 门 关 闭 时 间 对 分 义 管 末 
端 压 力 变 化 的 影响 较 大 。 
3.3 ”泄漏 位 置 的 影响 
为 了 探究 不 同 泄漏 位 置 对 分 又 管 阀门 末端 压力 
影响 , 取 下 游 分 又 管 的 长 度 Li = 六 =500 m.、 汇 漏 
量 为 0 的 10% ,保持 其 他 参数 不 变 ,数值 模拟 了 快 
速 关闭 阀门 和 绥 慢 关闭 阀门 两 种 方式 下 管道 的 泄漏 


3.3.1 快速 关闭 阀门 时 泄漏 位 置 的 影响 


采取 线性 快速 关闭 阀门 的 方式 ,阀门 的 关闭 时 
间 为 0.5 s, 对 不 同 的 泄漏 位 置 进行 数值 模拟 。 泄 漏 
位 署 分 别 位 于 直 管 段 距 上 游 水 库 200 .400 .600 .800 m 


从 数值 模拟 了 下 游 分 又 管 压力 变化 曲线 ,如 图 9 
所 未 。 
本 120 一 一 泄漏 点 距 上 游 200 m 
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(b) 阀门 末端 压力 变化 详 图 
到 9 不 同 泄漏 位 置 对 阀门 末端 压力 的 影响 


Fig.9 The influence of different leakage position 
on the pressure at the end of the valve 
由 图 9 (a) 可 以 看 出 ,在 不 发 生 泄漏 时 , 阀门 末 
端的 水 击 压力 峰值 最 大 , 随 着 泄漏 的 发 生 ,泄漏 位 置 
距离 上 游 越 远 ,阀门 末端 的 水 击 压力 峰值 越 小 ,但 在 
第 一 相 稳 定 区 域 ,不 同 的 泄漏 位 置 的 压力 峰值 大 致 
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相同 。 由 图 9(b) 可 以 得 到 ,在 距离 上 游 水 库 小 于 
600 m 位 置 处 泄漏 时 ,水 击 压力 曲线 呈现 凸 状 , 当 距 
离 上 游 水 库 大 于 800 m 处 发 生 泄 漏 时 ,压力 曲线 呈 
现 凹 状 。 因 此 ,可 以 根据 压力 曲线 呈现 的 凸凹 状 , 大 
致 判断 泄漏 所 发 生 的 位 置 ,减少 泄漏 检测 的 工作 量 ， 
为 进一步 的 精准 检测 定位 奠定 了 基础 。 
3.3.2 缓慢 关闭 阀门 时 泄漏 位 置 的 影响 

采取 线性 缓慢 关闭 阀门 的 方式 ,阀门 关闭 时 间 
为 4s, 对 不 同 泄漏 位 置 进 行 数值 模拟 。 汇 漏 位 置 位 
于 直 管 段 距离 上 游 水 库 100 .200 .300 .400 、500 .600、 
700 .800 .900 m 处 ,数值 模拟 了 阀门 末端 压力 变化 曲 
线 如 图 10 所 示 。 
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(b) 压力 曲线 面积 积分 
图 10 不 同 泄漏 位 置 对 阀门 末端 压力 的 影响 
Fig. 10 The influence of different leakage position 


on the pressure at the end of the valve 

由 图 10(a) 可知 ,水 击 波 压力 峰值 在 第 一 相 4 s 
时 达到 最 大 , 随 着 泄漏 位 置 距离 上 游 水 库 越 来 越 远 ， 
所 产生 的 水 击 压强 峰值 越 小 ,压力 尖 点 处 也 变 的 越 
来 越 平滑 ,水 击 波 衰减 速率 也 渐渐 放 缓 。 压 力 变 化 
曲线 横 、 纵 坐标 所 围 成 的 面积 可 以 得 到 压力 曲线 面积 
积分 ,由 图 10(b) 可知 ,泄漏 位 置 越 靠近 上 游 , 所 产生 
的 压力 曲线 面积 积分 越 大 ,在 距离 上 游 水 库 100 m 
时 ,达到 最 大 压力 曲线 面积 积分 为 877.943 m， s, 当 
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泄漏 位 置 距离 上 游 水 库 900 m 时 ,达到 最 小 压力 曲 
线 面积 积分 为 336.38 m . s, 并 且 随 着 时 间 的 推移 ， 
各 线 间 的 面积 积分 差 值 逐渐 增 大 。 因 此 , 当 管 道 发 
生 泄 漏 时 ,可 以 通过 压力 曲线 的 面积 积分 来 判断 泄 
漏 所 发 生 的 位 置 。 


3.4 ”泄漏 量 的 影响 


为 了 探究 不 同 泄漏 量 对 分 又 管 阀 门 末 端 压 力 的 
影响 , 取 下 游 分 又 管道 的 长 度 = /=500 m、\ 汇 漏 
位 置 距离 上 游 水 库 500 m 处 ,保持 其 他 参数 条 件 不 
变 ,数值 模拟 了 快速 关闭 阀门 和 缓慢 关闭 阀门 两 种 
方式 下 管道 的 泄漏 过 程 。 

3.4.1 快速 关闭 阀门 时 泄漏 量 的 影响 

二 采取 线性 快速 关闭 阀门 的 方式 ,阀门 的 关闭 时 
澡 0.5 s, 对 不 同 的 泄漏 量 进行 数值 模拟 , 泄漏 量 
放 


取 为 0.0. 1% Qo \1% Qo、10% Oo 、20% Q6o ,得 到 
了 EPB 游 分 又 管 压力 变化 曲线 如 图 11 所 示 。 
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图 11 不同 汇 漏 量 对 下 游 阀门 末端 压力 曲线 的 影响 


"= Fig.11 
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End pressure curve of downstream valve 


ch 


with different leakage 

由 图 11 可 知 , 随 着 泄漏 量 的 增 大 ,各 个 泄漏 量 
之 间 压 力 峰 值 的 差 值 就 越 大 ,在 第 3 个 周期 处 达到 
的 最 大 差 值 为 24.91 mm ,并 且 随 着 时 间 的 推移 ,衰减 
速率 在 逐渐 变 慢 , 当 汽 漏 量 为 0.0.1% 0 、1% Qo 时 ， 
由 于 泄漏 量 较 小 ,因此 压力 曲线 的 变化 并 不 明显 ,最 
大 压力 峰值 差 值 只 有 3.42 m, 但 当 泄漏 量 为 10% @u、 
20% Oo 时 压力 曲线 变化 差异 十 分 明显 ,但 不 同 的 泄 
漏 量 所 表现 出 的 衰减 趋势 较为 一 致 ,因此 ,在 较 小 泄 
漏 量 下 并 不 能 察觉 出 压力 变化 的 明显 差异 。 
3.4.2 缓慢 关闭 阀门 时 泄漏 量 的 影响 

采取 线性 缓慢 关闭 阀门 的 方式 ,阀门 的 关闭 时 
间 为 4 s, 泄漏 量 分 别 取 为 0.0. 1% 0 .0.5% 6、 
1% 0 53% Qo、10% 0 15% Qo、20% 0u ,数值 模拟 了 
下 游 分 义 管 压力 变化 曲线 ,如 图 12 所 示 。 
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12 不 同 泄漏 量 对 下 游 阀门 末端 压力 曲线 的 影响 


Fig.12 Pressure curve of downstream valve end 


器 


with different leakage 

由 图 12(a) 可知, 当 汇 漏 量 为 0\0.1% Oo .0.5% Qo 
时 ,阀门 末端 压力 信号 变化 差异 不 大 ,但 随 着 泄漏 量 
的 增 大 ,阀门 末端 压力 峰值 在 逐渐 减 小 ,衰减 速率 逐 
渐变 慢 ,在 4s 时 不 同 泄漏 量 的 水 击 压力 峰值 达到 最 
大 ,最 大 值 为 68.35 m, 各 线 间 压力 峰值 的 差 值 随时 
间 逐 渐 增 大 ,4 s 时 差 值 最 小 为 8. 50 m; 由 图 12(b) 
可 知 , 随 着 泄漏 量 逐 渐 增 大 ,阀门 末端 压力 曲线 面积 
积分 在 逐渐 减 小 ,并 且 面 积 积分 的 增长 速率 逐渐 减 
小 。 当 泄漏 量 为 0.1% Qo 时 ,达到 最 大 压力 曲线 面 
积 积分 为 973.81 ms, 当 泄漏 量 为 20% Oo 时 ,达到 
最 小 压力 曲线 面积 积分 为 504. 094 m. s, 并 且 随 着 
时 间 的 增长 ,各 个 泄漏 量 之 间 压 力 曲 线 面 积 积分 差 
值 越 来 越 大 。 


4 结 论 


在 复杂 地 形 的 输 水 系统 中 ,分 勾结 构 输 水 管道 
汇 漏 瞬 变 特性 引起 了 大 量 学 者 的 关注 。 本 研究 构建 
了 分 又 结构 输 水 管道 泄漏 数学 模型 ,并 验证 了 模型 
的 正确 性 ,所 得 结果 如 下 。 
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1) 不 同 分 又 管 长 度 对 阅 门 未 端 压力 峰值 的 影响 
较 大 ,分 叉 管 长 度 越 长 ,阀门 末端 的 压力 峰值 越 小 ， 
并 且 在 幅 值 处 的 变化 越 为 平缓 , 训 减 速率 也 逐渐 
变 慢 。 

2) 压 力 曲 线 的 差异 随 阀门 关闭 时 间 增加 而 增 
大 , 尖 点 处 压力 峰值 更 为 平滑 ,峰值 越 低 ,在 2。 闪 站 
关闭 时 间 内 ,压力 曲线 周期 的 变化 并 不 明显 。 

3 ) 随 着 泄漏 量 增 大 ,阀门 未 端的 压力 峰值 越 小 ， 
在 小 于 1% 0u 的 泄漏 量 下 并 不 会 引起 压力 信号 的 明 
变化 ,但 不 同 的 泄漏 量 所 引起 的 压力 变化 在 趋势 
上 较为 一 致 。 

4) 管 道 的 泄漏 位 置 对 水 击 有 较 大 影响 ,泄漏 位 
置 越 靠近 上 游 ,水 击 压力 峰值 越 大 ,并 且 压力 曲线 在 
幅 信 处 会 出 现 凸 起 趋势 ,但 当 泄漏 位 置 在 管道 下 半 
段 醒 ,压力 曲线 在 幅 值 处 会 出 现 止 陷 趋 势 ,因此 可 以 
柱 据 压力 曲线 变化 趋势 判断 泄漏 所 发 生 的 位 置 。 
马 分 又 结构 答 水 管道 泄漏 瞬 变 流 特性 的 系统 研 
EEN 有效 的 改进 了 管道 洪 漏 检测 方法 ,提高 了 泄漏 位 
锯 的 精度 ,为 保障 输 水 管道 的 安全 稳定 运行 提供 了 
有 定 的 理论 支撑 。 
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